n -/-<H 


NASA Army Research Laboratory 

Technical Memorandum 106623 Technical Report ARL-TR-475 

AIAA-94-2937 


Vibration Signature Analysis of a Faulted 
Gear Transmission System 


F.K. Choy and S. Huang 
University of Akron 
Akron, Ohio 

J.J. Zakrajsek 
Lewis Research Center 
Cleveland, Ohio 

R.F. Handschuh 
Vehicle Propulsion Directorate 
U.S. Army Research Laboratory 
Lewis Research Center 
Cleveland, Ohio 

D.R Townsend 
Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

30th Joint Propulsion Conference 

cosponsored by the ALAA, AS ME, S AE, and AS EE 

Indianapolis, Indiana, June 27-29, 1994 



National Aeronautics and 
Space Administration 


U.S. ARMY 



RESEARCH LABORATORY 


VIBRATION SIGNATURE ANALYSIS OF 
A FAULTED GEAR TRANSMISSION SYSTEM 


F. K. Choy and S. Huang 
Department of Mechanical Engineering 
University of Akron 
Akron, Ohio 44325 

J.J. Zakrajsek 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

R. F. Handschuh 
Propulsion Directorate 
U.S. Army Aviation Systems Command 
Lewis Research Center 
Cleveland, Ohio 44135 

and 

D.P. Townsend 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


Abstract 

A comprehensive procedure in predicting faults in gear 
transmission systems under normal operating conditions is 
presented. Experimental data was obtained from a spiral 
bevel gear fatigue test rig at NASA Lewis Research Center. 
Time synchronous averaged vibration data was recorded 
throughout the test as the fault progressed from a small 
single pit to severe pitting over several teeth, and finally 
tooth fracture. A numerical procedure based on the Wigner- 
Ville distribution was used to examine the time averaged 
vibration data. Results from the Wigner-Ville procedure are 
compared to results from a variety of signal analysis 
techniques which include time domain analysis methods and 
frequency analysis methods. Using photographs of the gear 
tooth at various stages of damage, the limitations and 
accuracy of the various techniques are compared and 
discussed. Conclusions are drawn from the comparison of 
the different approaches as well as the applicability of the 
Wigner-Ville method in predicting gear faults. 

I. Introduction 

In the last two decades, with the ever increasing 
demand for higher power and performance, premature 
failures in high performance machinery often result in major 
production and financial losses, and sometimes even lead to 
catastrophic consequences. In the aerospace industry, one of 
the major concerns is the rotorcraft transmission system. 
Currently, on-board operational condition monitoring 


systems for rotorcraft transmissions often fail to provide 
sufficient time between warning and failure so that safety 
procedures can be implemented. In addition, inaccurate 
interpretation of operating conditions often result in false 
alarms and unnecessary repairs and downtime. The develop- 
ment of an accurate machine health monitoring system for 
rotorcraft transmission systems will increase safety and 
reduce financial investments. 

Three major types of vibration analysis have been used 
in the early detection of machine failure, namely, frequency 
analysis, various techniques using the time domain analysis, 
and joint time-frequency analysis. 

For frequency analysis, the amplitude spectrum of the 
measured vibration signal is calculated and displayed in a 
continuous manner. 15,18 This method is particularly power- 
ful in identifying vibrations due to various mechanical 
elements in the system at their operating frequencies. For 
example, a gearbox will experience vibration at the funda- 
mental and multiple harmonic frequencies of the gear tooth 
meshing frequency, which are usually different for each pair 
of gears in the gearbox. The changes in the frequency 
spectrum, such as an increase in the vibration amplitude at 
a particular frequency, or the appearance of modulation 
sidebands about that frequency, can often be identified and 
related to a single element in the complex system. Further- 
more, abnormal vibration signals, caused by wear, deforma- 
tion of the machine element, "ghost component” etc., which 
are sometimes difficult to quantify in the time domain, can 
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be readily identified in the frequency spectrum. However, 
the spectral analysis technique is difficult to apply in a 
highly complex system, such as a helicopter gearbox, which 
may have as many as 30 gears and 50 bearings. The large 
number of spectral lines which these elements produce 
during normal operation can make it difficult to even detect 
changes in the spectrum, let alone to identify the sources by 
which these changes are generated. In addition, the modula- 
tion of the vibration signal of early fatigue cracks in the 
mechanical elements often yield sidebands which extended 
over a broad frequency range with low amplitude. Under 
these circumstances, the background noise level in the 
spectrum can easily mask the presence of the sidebands, and 
mechanical defects in the system may go undetected. 

Time domain analysis is also used for the early 
detection of machine faults. 9-11,17 ’ 19 " 21 The digitized time 
domain vibration signal is averaged over a large number of 
cycles, synchronous with the running speed of a particular 
machine element. This time synchronous averaging removes 
the background noise, leaving a more accurate estimate of 
the vibration signal of the component of interest. This 
process can be repeated for each machine element in the 
system to provide the vibration signal relative to that 
particular element. A variety of statistical based techniques 
are then applied to the averaged signal to extract specific 
fault information. In most machinery, vibration signals are 
often characterized by a variation of both amplitude and 
frequency with respect to time. For most mechanical faults, 
both the amplitude and frequency contents of the signal will 
shift during the period in which the damage is generated. 
The nature of these short-term variations in the time signal 
could prevent an accurate indication of the characteristics of 
the mechanical fault present. 

The joint time-frequency analysis approach applies the 
Wigner-Ville Distribution (WVD) 1_3 ’ 14 ’ 16 on the time 
vibration signal of the system. Unlike the Fourier Transform 
process, the WVD provides a instantaneous frequency 
spectrum of the system. The spectral density of the funda- 
mental exciting frequency and its sidebands changes as the 
shaft rotates through a complete revolution. If the vibration 
repeats itself each revolution, a Fourier Transform will show 
a constant frequency spectrum, while the changing spectral 
density with time from the WVD provides information 
concerning the frequency distribution concentrated at that 
instant around the instantaneous frequency. Some success 
has been achieved in applying the WVD in gear transmis- 
sion systems 5,6,8,12,16 to recognize faults at various locations 
of the gear. 


The time domain analysis includes the time synchronous 
averaged signal, and the statistical based techniques FMO, 
FM4, NA4*, and NB4* applied to the time averaged signal. 
The joint time-frequency analysis uses the WVD on the 
vibration data with special windowing techniques. All of the 
analysis methods are applied to the vibration data of a 
failure from the spiral bevel gear fatigue test rig in the 
NASA Lewis Research Center. Results from the various 
methods are compared and general conclusions are drawn 
from the results. 

II. Technical Approach 

As discussed in the previous section, three major 
methodologies: A) the frequency domain approach, B) the 
time domain approach, and C) the joint time-frequency 
approach are used in this study. The following is a descrip- 
tion of the three methodologies: 

(A) Frequency Domain Techniques 

The frequency spectrum is found by applying a discrete 
FFT on the time averaged signal x(t), such that the spectral 
components are 

X(k) = T l x(t) expLi^i (!) 

to N 


x(t) = time domain signal 

X(k) = frequency domain signal 

T = sampling time interval 

N = number of data points 

The frequency components are examined in the frequency 
domain and compared with those obtained at various stages 
of the fault development in the spiral bevel pinion. 

(B) Time Domain Techniques 

Four different time domain techniques for early 
detection of gear tooth damage are used in this study for 
evaluation and comparison. All of the time domain tech- 
niques are applied to the vibration signal after it has been 
time synchronously averaged. These techniques are: FMO, 

FM4, NA4*, and NB4*. These parameters are defined as 
follows: 

FMO 17 is a course fault detection parameter that 
compares the maximum peak to peak amplitude to the sum e 
of the mesh frequencies and its harmonics 


where 


The objective of this paper is to examine and compare 
the three different approaches to detect gear wear and 
failure. The frequency domain analysis is based on the 
spectral display from a Fast Fourier Transform algorithm. 
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where 


PP = maximum peak to peak amplitude in signal 
A(fj) = amplitude at mesh frequency and harmonics 
n = total number of harmonics in frequency range 


FM4 17 is an isolated fault detection parameter, and is 
given by the normalized kurtosis, of the resulting difference 
signal as 


FM4 = 



(dj - d) 4 


r N -,i2 

l (di - d) 2 


(3) 


where 

d(t) = A(t) - R(t) 

d = mean value of d(t) 

A(t) = original time synchronous signal 
R(t) = regular meshing components plus their first 
order side bands 

N = total number of data points in the time signal 


NA4* 4,20,21 is a general fault detection parameter with 
trending capabilities. A residual signal is constructed by 
removing regular meshing components from the time 
averaged signal. For NA4, the first order sidebands stay in 
the residual signal and the fourth statistical moment of the 
residual signal is then divided by the averaged variance of 
the residual signal, raised to the second power. The average 
variance is the mean value of the variance of all previous 
records in the run ensemble. This allows NA4 to compare 
the current gear vibration with the baseline of the system 
under nominal conditions. NA4 is given by the quasi- 
normalized kurtosis equation shown below: 


NA4*(M) = 


n£ (r, - r) 4 
i=l 



(4) 


where 

r = residual signal 

r = mean value of residual signal 

N = total number of data points in one time record 

i = data point number in time record 
j = time record number 

M = current time record number in run ensemble 


This provides NA4 with enhanced trending capabilities, in 
which the kurtosis of the current signal is compared to the 
variance of the locked baseline signal under nominal 
conditions. 


NB4 is another parameter similar to NA4 that also uses 
the quasi-normalized kurtosis given in equation (4). The 
major difference is that instead of using a residual signal, 
NB4 uses the envelope of the signal banpassed about the 
mesh frequency. Again, as with NA4*, NB4* is an enhance- 
ment to the NB4 parameter, in which the value of the 
average variance is "locked” when the instantaneous 
variance exceeds a pre-determined value. The equation for 
NB4* is given below: 


NB4*(M) = 


N 

n£ (sj - I) 4 

i=l 


fl N 

r N 
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i < s y - *j) 2 


l M j=l 

Li=i J 
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and 

s(t) = magnitude of {b(t) + i {H[b(t)]}) 


(5) 


( 6 ) 


where 

b(t) 

H[b(t)] 

N 


j 

M 


= time averaged signal bandpassed filtered 
about the meshing frequency 
= the Hilbert Transform of b(t) 

= total number of data points in one time 
record 

= data point number in time record 
= time record number 

= current time record number in run ensemble 


(C) Joint Time-Frequency Technique 

To examine the vibration signal in a joint time- 
frequency domain, the Wigner-Ville method 1-3,14 is used in 
this study. While the FFT technique (eq. 1) can provide the 
spectral contents of the time signal, it cannot distinguish 
time phase change during a complete cycle of operation. In 
other words, it assumes that the time signals are repeatable 
for each time data acquisition window without considering 
the effects of any magnitude and phase changes during the 
sampling period. The Wigner-Ville distribution will provide 
an interactive relationship between time and frequency 
during the period of the time data window. The WVD 
(Wigner-Ville Distribution) can be written as: 

W(t, f) = £ xjt 4 i.1 x*jt - .1} e‘ j27rfx dt ( 7 ) 

• eo 2 2 


or in a discrete form as 
L 

W x (nT, f) = 2T £ x(nT* iT) x*(nT- iT) . e 

i=-L 


An enhancement to this parameter is given by NA4*, in 
which the value of the averaged variance is "locked" when 
the instantaneous variance exceeds a pre-determined value. 4 
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( 8 ) 


where 

W(t, f) = the Wigner-Ville distribution in both the 
time domain t and frequency domain f. 
x(t) = the time signal 

T = the sampling interval 

L = the length of time data used in the transform 

To allow sampling at the Nyquist rate and eliminate the 
concentration of energy around the frequency origin due to 
the cross product between negative and positive frequen- 
cy, 1 3 the analytic signal was used in evaluating the WVD. 
The analytic signal s(t) is defined as 

s(t) = x(t) ♦ jH[x(t)l (9) 

where 

H[x(t)] = the Hilbert transform of x(t) defined by: 



However, an alternative approach can be used to calculate 
the analytic signal using the frequency domain definition. 
The analytic signal s(t) can be evaluated by calculating the 
FFT of the time signal x(t), then setting the negative 
frequency spectrum to zero. The analytic signal can be 
obtained by evaluating the inverse FFT of the spectrum. To 
simplify the computational effort, the WVD can be evaluat- 
ed using a standard FFT algorithm. Adopting the convention 
that the sampling period is normalized to unity, equation (8) 
can be rewritten as 

L 

W x (n, f) = 2 £ x(n* i) x*(n- i) . e‘ j47tfl ( n > 
i=- L 

As for the continuous time case, it is necessary only to 
evaluate the WVD at time zero. 

Hence 

L 

w x (0, f) = 2 Y k(i) e ' j47rfi ( 12 > 

i=- L 

where 

k(i) = s(i) s* (-i). 

Equation 12 can be evaluated using the discrete FFT 
algorithm. 

In order to avoid a repetition in the time domain WVD, 
a weighting function 12 is added to the time data before the 


evaluation process. Such a process may decrease the 
resolution of the distribution, but it will eliminate the 
repetition of peaks in the time domain and the interpretation 
of the result will be substantially easier. 

m. Description of Experimental Procedure 

The fatigue damage on the test pinion shown in Figures 
1 to 8 was obtained using the spiral bevel gear fatigue test 
rig illustrated in figure 9, at the NASA Lewis Research 
Center. The primary purpose of this rig is to study the 
effects of gear tooth design, gear materials, and lubrication 
types on the fatigue strength of aircraft quality gears. 7 
Because spiral bevel gears are used extensively in helicopter 
transmissions to transfer power between nonparallel inter- 
secting shafts, the use of this fatigue rig for diagnostic 
studies is extremely practical. Vibration data from an 
accelerometer mounted on the pinion shaft bearing housing 
was captured using a personal computer with an analog to 
digital conversion board and anti-aliasing filter. The 
12-tooth test pinion, and the 36-tooth gear have:, 0.5141 in 
pitch, 35 degree spiral angle, 1 in. face width, 90 degree 
shaft angle, and 22.5 degree pressure angle. The pinion 
transmits 720 hp at nominal speed of 14,400 rpm. The test 
rig was started and stopped several times for gear damage 
inspection. The test was ended at 17.79 operational hours 
when a broken portion of a tooth was found visually during 
one of the shutdowns. 

IV. Discussions of Results 

A series of pictures showing the deterioration of the 
pinion teeth at various stages of the test are illustrated in 
figures 1 to 8. In figure 1 the initiation of a small pit on one 
of the pinion teeth during the first shutdown, at about five 
and a half hour into the test is shown. As the test pro- 
gressed, the rig was shut down seven more times to exam- 
ine the severity of the pitting and its relationship with the 
corresponding vibrations. Figures 2 to 4 show the increase 
of the damaged area at the pinion tooth as the elapsed time 
increased to 6.55, 8.55, and 10.03 hr respectively. Note that 
in figure 5, as the elapsed time increased to 12 hr, the 
damage of the pinion tooth increased to 75 percent of the 
tooth surface. At this stage, pitting also initiated on the 
adjacent tooth and continued to grow 21 s the time increased 
to 14.53 hr, figure 6. At 16.16 hr, the damage has grown to 
three adjacent teeth as shown in figure 7. The test was 
terminated when a breakage is detected on one of the three 
heavily pitted teeth at 17.79 hr, as shown in figure 8. 

Figure 10 depicts the running speed of the test rig 
during various stages of the experiment. Note that there is 
some fluctuations present in the running speed after each 
shutdown, with a magnitude of approximately 6 percent 
about the nominal pinion speed of 14,400 rpm. There is a 
sharp change in speed at approximately 8.75 hr. These 
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variations in speed create a substantial effect on the vibra- 
tion signal, which is amplified in the NA4 (fig. 13), NB4 
(fig. 14), and the WVD (fig. 21) analysis. 

Figure 1 1 shows the results of the FMO analysis. As 
seen from the figure, FMO shows only moderate changes as 
the damage starts and progresses. It does not provide any 
indication as the damages extends to the adjacent teeth, 
resulting in pinion tooth fracture. The majority of the 
variations in the FMO parameter are most probably due to 
the speed changes experienced during the test. 

Results from the FM4 parameter, as seen in figure 12, 
shows a possible reaction as the pitting started to occur, 
however, it does not provide any coherent indication of the 
severity of the pitting as the damage increased. In addition, 
it does not provide any information to distinguish the pitting 
of a single or multiple teeth. 

Results from NA4 and NA4* are illustrated in fig- 
ure 13. It is obvious from the figure that both NA4 and 
NA4* provide a very good indication of the pitting develop- 
ment on the pinion tooth. The magnitude of the parameter 
increases to a nondimensional value of 7 after shutdown #2 
at 6.55 hr, and further to a value of 17 when the pitting 
covers 75 percent of the tooth surface at 8.5 hr. As expect- 
ed, the ’’locked” denominator in NA4* provides a more 
robust indication as the pitting progresses. 4,21 Again, both 
parameters are very sensitive to the speed variations, 
especially after shutdown #3 and #6. 

The NB4 and NB4* parameters, as shown in figure 14, 
show a very similar trend to those of NA4 and NA4*, with 
a more robust indication to the severity of the damage. 
However, both NA4 and NB4 did not provide any type of 
indication as the damage spread to other teeth, and finally 
as tooth fracture occurred. 

Figure 15 shows a WVD (and corresponding intensity 
scale) and below the WVD the frequency spectrum (from 
FFT analysis) of a uniform sine wave signal shown at the 
left side of the figure. Note that only the frequency com- 
ponent of the input frequency is detected by both the 
frequency domain analysis and the joint time-frequency 
distribution. The WVD does not exhibit any changes during 
the 1 cycle (0 to 360 degrees) rotation of the shaft. When 
a short term amplitude and phase change is added to the 
system, as shown in figure 16, the frequency spectrum 
remains virtually unchanged. The WVD shows a dramatic 
change of the energy distribution pattern at the location 
where the change occurs. The lighter shades of the distribu- 
tion display indicates a smaller vibration amplitude, which 
is shown by the time signal at the left side of the figure. 
Such an effect could be possibly caused by a chipped or 
cracked tooth. Figure 17 shows the effects of a short term 
amplitude increase in the time signal to simulate vibrations 
caused by gear tooth surface damage. Note that the frequen- 


cy spectrum remains the same showing only the component 
of the exciting frequency while the WVD again provides a 
good indication of the amplitude increase by the widening 
of the shaded area to a diamond shape at the corresponding 
’’damaged” tooth location. Figure 18 shows the effects of a 
time decaying short term amplitude and phase change 
signal. The WVD shows a half diamond shape of shaded 
area, similar to that of figure 16, at the location where the 
amplitude and phase changes are presented. As seen in 
figure 18, the frequency spectrum gives very litde indication 
of the signal change. 

Figures 19 to 26 show the WVD and frequency spectra 
of the spiral bevel pinion vibration at various stages of 
damage, corresponding to the photographs given in figures 1 
to 8. Figure 19 shows the occurrence of the initial pitting at 
around 200 degrees from the triggering point of the gear, at 
5.5 hr, of running time. At 6.55 hr, the pitting on the tooth 
surface progressed to a more noticable stage, as shown in 
figure 2, the WVD pattern, figure 20, begins to adopt those 
of a short term amplitude and phase change as illustrated in 
figure 18. Due to the speed increase at 8.55 hr, the overall 
WVD amplitude increases substantially as shown in figure 
21. At the running time of 10.03 hr, the corresponding 
WVD in figure 22 shows the initiation of a cross pattern as 
surface pitting in the damaged tooth becomes more pro- 
nounced which can also be evidented from the photograph 
of the pinion gear in figure 4. This phenomenon is also 
evident in the frequency spectrum with the existence of 
sideband components. At 12.03 hr, when the pitting on the 
pinion tooth hax extended to about 75 percent of the tooth 
surface as shown in figure 5, the WVD pattern exhibits a 
solid cross pattern extending over the mesh freqeuncy and 
several of its adjacent sidebands. The high concentration in 
the WVD energy and the initiation of a second cross pattern 
at 14.53 hr, figure 24, shows the advancement of the pitting 
process on second tooth, as illustrated in figure 6. This is 
further confirmed by the large amplitude of sideband 
component(above mesh) in the frequency spectrum. At 
16.16 hr, as seen in figure 25, the WVD pattern changes, 
showing more advanced damage pattern similar to the 
multiples of the decay of a single short term amplitude 
increase and phase change demonstrated in figure 18. Such 
phenomenon is due to the pitting of three consecutive teeth 
in the pinion gear as shown in the picture given in figure 7. 
The frequency spectrum in figure 25 shows a substantial 
amplitudes increase in the sideband components. The 
discontinuity of the WVD at the mesh frequency, shown in 
figure 26, similar to the example shown in figure 16 due to 
the short term phase change, is probably the result of the 
instantaneous phase change caused by the fractured tooth, as 
illustrated in figure 8. The two cross patterns in the WVD 
is very distinct as the effects of pitting at the teeth adjacent 
to the fracture tooth become more pronounced. Note, also 
that, as given in figure 26, the amplitude of the sideband 
frequencies(above and below the mesh frequency) increase 
substantially. 
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V. Conclsuions 


6 . 


Based on the results of the application of the various 

aforementioned methods, the following conclusions can be 

made: 

1) The FMO parameter shows only moderate 
changes as the damage starts and progresses. 
It also fails to indicate the fracture of the 
pinion tooth. 

2) The FM4 parameter shows a possible reaction 
to the start of the pitting process, however, no 
coherent indication is provided for the growth 
and severity of the pit 

3) The NA4* and NB4* parameters show good 
reactions to the initial pitting damage and 
very nice indications for the growth and 
severity of the pitting damage. However their 
indications for the tooth fracture is somewhat 
unclear. 

4) The WVD provides vital information con- 
cerning both the severity and the location of 
the pitting process in the gear system. 

5) The fracture of the gear tooth and its exact 
location can be pinpointed using the WVD 
technique. However a machine vibration 
signature database is required to interpret the 
resulting WVD. 

6) The occurrence and the severity of gear tooth 
failure can be reliably detected using a com- 
bination of the time averaging, the frequency 
analysis, and the WVD techniques. 
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Figure 1 . — Photograph of pinion damage at t = 5.50 hr. 


Figure 2. — Photograph of pinion damage at t = 6.55 hr. 




Figure 3. — Photograph of pinion damage at t = 8.55 hr. 
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Figure 5.— Photograph of pinion damage at t = 12.03 hr. Figure 6.— Photograph of pinion damage at t = 14.53 hr. 



Figure 7. — Photograph of pinion damage at t = 16.16 hr. Figure 8. — Photograph of pinion damage at t = 17.79 hr (end of run). 
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Figure 9. — Spiral bevel gear fatigue test rig at NASA Lewis Research Center. 
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Figure 1 1 FMO parameter at various run time. 
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Figure 12.— FM4 parameter at various run time. 
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Figure 13.— NA4 and NA4* at various run time. 



Figure 15.— Example of WVD on a regular sine wave time signal. 
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Figure 16.— Example of WVD on short term amplitude and phase change in a signal. 
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Figure 17. — Example of WVD on short term amplitude increase in a signal. 



Figure 18.— Example of WVD on time decay of a short term amplitude and phase change in 
a signal. 
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Figure 19.— Application of WVD to the spiral bevel gear vibration signal at 5.50 hr. 
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Figure 21 .—Application of WVD to the spiral bevel gear vibration signal at 8.55 hr. 



Figure 20. — Application of WVD to the spiral bevel gear vibration signal at 6.55 hr. 



Figure 22. — Application of WVD to the spiral bevel gear vibration signal at 10.03 hr. 







Figure 23.— Application of WVD to the spiral bevel gear vibration signal at 12.03 hr. 



Figure 25.— Application of WVD to the spiral bevel gear vibration signal at 1 6.1 6 hr. 



Figure 24.— Application of WVD to the spiral bevel gear vibration signal at 14.53 hr. 



Figure 26.— Application of WVD to the spiral bevel gear vibration signal at 1 7.79 hr. 
(end of run). 
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